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It is well known that the Casimir energy of bulk fields induces a non-trivial potential for 
the compactification radius of higher- dimensional field theories. On dimensional grounds, 
the 1-loop potential is ~ l/R 4 . Since the 5d gauge coupling constant g 2 has the dimen- 
sion of length, the two-loop correction is ~ g 2 /R 5 . The interplay of these two terms leads, 
under very general circumstances (including other interacting theories and more compact 
dimensions), to a stabilization at finite radius. Perturbative control or, equivalently, a 
parametrically large compact radius is ensured if the 1-loop coefficient is small because 
of an approximate fermion-boson cancellation. This is similar to the perturbativity ar- 
gument underlying the Banks-Zaks fixed point proposal. Our analysis includes a scalar 
toy model, 5d Yang- Mills theory with charged matter, the examination of S l and S l /Z% 
geometries, as well as a brief discussion of the supersymmetric case with Scherk-Schwarz 
SUSY breaking. 2-Loop calculability in the S l / Z2 case relies on the log-enhancement of 
boundary kinetic terms at the 1-loop level. 



1 Introduction 



Higher-dimensional field theories arise in the low energy limit of string- or M-theory, 
which is our best candidate for a theory of quantum gravity. Independently, compacti- 
fied higher-dimensional models provide many interesting possibilities for the unification 
of known fields and interactions. Familiar examples are the appearance of 4-dimensional 
gauge theories as a manifestation of higher-dimensional diffeomorphism invariance or the 
unification of known bosons and fermions in higher-dimensional multiplets of supersym- 
metry (SUSY). Thus, we consider higher- dimensional compactified models a promising 
ingredient in possible physics beyond the standard model, which makes the further inves- 
tigation of stabilization mechanisms for the compactification radius an interesting and 
potentially important subject. 

One very generic ingredient in the dynamics of the compactification radius is the 
Casimir energy of massless bulk fields [1]. As a simple example, consider 5d general 
relativity with a vanishing cosmological constant and a set of massless 5d fermions and 
bosons. Compactifying to 4d on an S 1 with physical volume 2irR, one finds a flat 4d 
effective potential for the radion R. This flatness is lifted by the 1-loop Casimir energy 
which, on dimensional grounds, is ~ 1/R 4 and does not lead to a stable finite- R solution. 1 
Clearly, to overcome the 1/R 4 behaviour which is too simple for stabilization, one has 
to introduce a mass scale into the potential, which can be achieved, for example, by 
considering warped compactifications [2], massive bulk matter or brane localized kinetic 
terms for bulk fields [3]. 

In this paper, we point out that the required mass scale is, in fact, generically present 
in the simplest interacting higher-dimensional field theories, such as \(f) 4 theories, gauge 
theories with coupling constant g, or models with Yukawa interactions. In 5d the above 
coupling constants are dimensionful, leading to a 2-loop Casimir energy contribution 
~ X/R 5 or g 2 /R 5 . Thus, radius stabilization will generically arise at the 2-loop level by a 
balancing of the 1/R 4 and the 1/R 5 contributions without the need to invoke any extra 
effects or operators. 2 Generically, the compactification scale is set by the lowest of the 
strong interaction scales of various 5d field theories present in a given model. 

We note that a different 2-loop stabilization mechanism was previously considered 
in the context of 6d A</> 3 theory, where the coupling is dimensionless and a logarithmic 
-R-dependence arises at the 2-loop level [5]. Furthermore, the possibility of 2-loop stabi- 
lization based on the vanishing of the 1-loop contribution 1/R 4 and a balancing of the 
1/R 5 and the \n(R)/R 5 terms has been pointed out in [6]. Two-loop corrections to the 
4d Casimir effect have been considered by many authors (see, e.g., [7]). 

The paper is organized as follows. In Sect. 121 the above idea is illustrated using the 
simple example of 5d \<fi 4 theory. We emphasize in particular that, by a judicious choice 

1 Introducing a non-zero 5d cosmological constant A§, a stable solution can be found by balancing the 
resulting 27ri?As contribution against the Casimir energy. However, a positive 4d cosmological constant 
results whose scale is set by R and which is therefore generically too large. 

2 Note that this differs qualitatively from the results of the early discussion of higher-loop Casimir 
stabilization in [4] (see below for more details). 
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of the field content, stabilization at moderately large radii, R 3> A, can be achieved, 
such that higher-loop corrections are negligible. This is similar to the way in which 
a perturbatively controlled non-trivial fixed point arises in the proposal of Banks and 
Zaks [8]. 

Section |H] extends the analysis to a 5d Yang-Mills theory with charged bosons and 
fermions. Amusingly, all 2-loop integrals reduce straightforwardly to the simple scalar 
case. Controlled 2-loop stabilization at large radius can be achieved, e.g., in the large N 
limit of SU(iV) gauge theory with appropriate matter content. Furthermore, it is shown 
that our 2-loop stabilization mechanism extends straightforwardly to SUSY models with 
Scherk-Schwarz SUSY breaking. 

In Sect. HI we provide a qualitative discussion of the phenomenologically more inter- 
esting cases of S l jZi and of the S 1 with 3-branes. Since the finite and calculable 2-loop 
bulk contribution mixes with the 1-loop effect induced by brane operators with unknown 
coefficients, a complete predictivity just on the basis of the field content can not be 
achieved. In the generic case, the 2-loop effect considered here represents an 0(1) correc- 
tion to the previously discussed stabilization by brane-kinetic terms. We point out the 
interesting and natural limit of logarithmically enhanced brane-localized gauge-kinetic 
terms, which allows one to neglect the bulk-2-loop effect and to achieve full predictivity 
just on the basis of the particle spectrum. 

A summary of our results as well as a discussion of possible further research direc- 
tions, in particular the applicability to SUSY models on S 1 /Z 2 and to the case of more 
than 5 dimensions, can be found in Sect. |SJ 

The evaluation of the relevant loop integrals and non-trivial SUSY-based checks of 
the 2-loop gauge theory calculation are described in two appendices. 



Consider the 5d theory of Einstein gravity and massless real scalar with classical action 



compactified on an S 1 with radius R. Clearly, a variation of the metric background field 
<?55 is equivalently described by a variation of the volume 2nR. In the following, we will 
always use a 5d Minkowski metric treating R as our volume or radion degree of freedom. 
It will not be necessary to perform a Weyl rescaling of the metric to manifestly separate 
graviton and radion degrees of freedom in the 4d effective theory. 

Assuming that Mp.5 ^> 1/A, we can consistently neglect gravitational interactions. 
On dimensional grounds, the effective potential for R then reads 



2 A example 
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where the are n-loop coefficients. Note that, even in the limit of large 5d Planck 
mass, has to include the 5d graviton contribution. 

Radius stabilization can be achieved already at the 2-loop level. Indeed, if is 
negative and positive, the 2-loop potential is minimized at 

However, the 4d cosmological constant at the minimum is negative. This can be reme- 
died either by adding a 3-brane with appropriately tuned positive tension or a 5d bulk 
cosmological constant. In the second case, an extra independent contribution to the 4d 
effective potential results, 

V CC (R) = 2ttA 5 R. (4) 

Requiring both V'(R) and V(R) to vanish at the same point determines the precise value 
of A5 and gives rise to a slightly shifted minimum at 3 



Unfortunately, assuming that = 0(1), it is immediately clear that higher-loop 
terms can not be neglected in the vicinity of the above 2-loop minimum. 4 This situation 
may, in fact, be generic. In this case we can not do more than express the justified hope 
that, in many models, higher-loop effects will significantly change but not destroy the 
minimum found at the 2-loop level. Our first conclusion is therefore that radius stabi- 
lization by higher-loop Casimir energy is presumably generic (in the sense of occuring in 
a large fraction of models) and the resulting compactification scale is of the order of the 
strong interaction scale of the most strongly coupled of the bulk field theories. 

However, in specific examples a quantitatively controlled minimum based on the 
2-loop approximation can occur. Indeed, in models where \c^\ <C |c^|, the 2-loop min- 
imum is at R ^> A and higher-loop effects are suppressed (assuming that no undue 
enhancement of the coefficients with n > 3 occurs). As a concrete realization, con- 
sider the 0(N S ) symmetric generalization of the above scalar A0 4 lagrangian, 

together with Nf fermions which are not (or only very weakly) coupled to the scalars. It 
is clear that, in this case, 

C W ~ 4jV7 _ N s - 5 , (7) 

3 Strictly speaking, since we are not working in the Einstein frame, the equation of motion for R is 
V + -KMp 5 1Z/t — 0, where 72-4 is the 4d curvature. Demanding V — at the minimum also ensures that 
minimization of V is equivalent to solving the equation of motion. 

4 We could, of course, improve our estimates by extracting appropriate loop suppression factors from 
the coefficients using naive dimensional analysis (see, e.g., [9]). However, since this would not affect 
our argument qualitatively, we do not enter such a more detailed discussion. Alternatively, one can 
imagine that these factors have already been absorbed in a redefined coupling. 
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(which is the difference of fermionic and bosonic degrees of freedom, including the 5 
graviton polarizations) while 

c (n) ~ iV s n for n > 2 . (8) 

Thus, taking N s large while keeping 4Nf—N s —5 ~ 0(1) and negative, one has R ~ XN^ 
at the 2-loop minimum and all higher-loop terms are suppressed. 

Finally, we now want to fill in the explicit numbers for the first two loop coeffi- 
cients used above. As already mentioned, is proportional to the difference of on-shell 
fermionic and bosonic degrees of freedom of the 5d theory, 

^ (-Afcrmions -Abosons) Cq , (9) 



where [1,3] (see also Appendix A) 

c« = (10) 

The 2-loop coefficient for a single scalar is due to the "figure-8" diagram, which has 
previously been derived using the winding mode expansion [6]. In our context, it is 
crucial that the tree-level masslessness in 5d is maintained by an appropriate 1-loop 
counterterm. 5 The result reads (for explicit calculations see Appendix A) 

c (2) = C(3) 2 , . 

8(2vr)9 ■ 1 j 

Going from a single scalar to the 0(A^)-symmetric model, the coefficient of Eq. (ITTT) 
has to be multiplied by + 2N s )/3. For our simple scalar example to work, it is 
important that c*- 2 -* > 0. 

Thus, a quantitatively controlled 2-loop minimum arises from the potential 

v ^ = WfR> { 3 « 5 >< 4iv ' ~ N > ~ 5 > + l§f^ + 2 ^ s} • < 12 > 

in the specific large- N s limit described above. 

While this simple analysis shows that it is quite easy to achieve stabilization in 
a scalar <f) A theory, a more interesting and realistic case is that of a 5d gauge theory. 
In particular, the masslessness assumption introduced above, which is unnatural for an 
interacting scalar, will be natural for gauge bosons and charged fermions. 



3 Gauge theory 

We now turn to the case of a 5d gauge theory compactified on S 1 with the action 

S = J d 5 x^ (^M 3 P , 5 n 5 - ±tT f {F MN F MN ) + ( J D M $)t( J D M $) + mp . (13) 

5 There is however a finite (nonlocal) positive correction to the scalar mass squared of the 4d zero 
mode. 
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Figure 1: Two-loop diagrams from the gauge sector. 





Figure 2: Two-loop diagrams from the matter sector. 

Here the bosonic and fermionic matter fields transform in some representation r of the 
gauge group G, the field strength is Fmn = — i[-Djw> D^} with Dm = 9m +iAm, and the 
trace is in the fundamental representation with generators normalized by 2trj(T a T fe ) = 
5 ab . 

The 2-loop vacuum energy contributions come from the diagrams depicted in Figs. [2 
and 121 It turns out that, by simple algebraic manipulations, they can all be reduced to 
the 2-loop integral encountered in the scalar model above. Therefore, the complete result 
can be expressed in terms of the constant 

c?) = H dim(G) (i4) 

where dim(G) is the dimension of the gauge group (the 1/8 of Eq. is a symmetry 
factor characteristic of the A0 4 model). 

Including the symmetry factors and numerator algebra of the various diagrams and 
accounting for the trace normalization through the constant T(r), defined by tr r (T a T 6 ) = 
T{r)5 ab (which is also known as the Dynkin index of the representation r), one finds 

cSL/oS = \d{d-l)T{a)-\{d-l)T{a) + \T{a) = +|T(a), 



4 2 LM 2) = dT{r)-\T{r) = +|T(r), (15) 

41ionM 2) = (2-d)T(r) = -3T(r). 
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The sum of these coefficients defines of Eq. (J2J) (with A replaced by g 2 ) and thereby the 
2-loop contribution to the Casimir energy arising from gauge fields and gauged matter. 
To facilitate comparison with other calculations, we have made explicit the ^-dependence 
before setting d = 5 and specified the contributions arising from the separate Feynman 

(2) 

gauge diagrams. Specifically, the three contributions to c^, ector come from the "figure 8" , 
the "setting sun" diagram, and the "setting sun" diagram with ghosts (in this order, 
cf. Fig. [T]). Analogously, the two contributions to c scaUr arise from the "figure 8" and 
the "setting sun" diagram (cf. Fig. |2J. A non-trivial check of these results based on 
the vanishing of the Casimir energy in models with unbroken SUSY can be found in 
Appendix B. 

One can see that in a pure gauge theory, one again encounters precisely the situation 
outlined above: S- 1 ' is negative while C ( 2 ) is positive so that 2-loop stabilization is au- 
tomatic. In order to maintain perturbativity, we need to reduce the one-loop coefficient 
without affecting the other c^ n \ This is most easily achieved by considering large groups 
(e.g. SU(iV) with N large, where 

c (n) _ 

iV 1+n ) and adding fermions uncharged under 
the gauge group to reduce to 0(1). This will lead to R ~ g 2 N 3 at the 2-loop min- 
imum and result in a relative suppression of the n-loop (n > 2) contribution near this 
minimum by 1/N 2n ~ 4 . In fact, this situation may arise fairly naturally since higher-loop 
coefficients are dominated by the most strongly coupled gauge group factor, so that it is 
sufficient to require the relevant fermions to be neutral only under this part of the group. 
We have thus seen that it is easy to stabilize the S 1 radius R in a controlled fashion. Of 
course, to build a realistic model one would have to include branes or to consider more 
sophisticated geometries allowing for chiral fermions. 

We close this section by commenting on a possible SUSY version of this scenario. 
Consider a model containing a 5d supergravity multiplet, Ny = dim(G) vector multiplets 
defining a super Yang-Mills theory with gauge group G, and Nh hypermultiplets in a 
representation r of G. Let us break SUSY from M = 2 to M = by introducing Scherk- 
Schwarz boundary conditions on the S 1 [10]. The effect of this breaking is a mass shift for 
gauginos, gravitinos and hyperscalars. Their Kaluza-Klein (KK) masses become m n R = 
n + u, where the real parameter u is known as the Scherk-Schwarz parameter. Already 
at this point it is clear that our stabilization mechanism is qualitatively unchanged: The 
Scherk-Schwarz parameter is dimensionless and our basic formula, Eq. (J2J), remains valid. 
Of course, the coefficients cW are now functions of uj (which vanish for lo — 0). 

The 1-loop contribution to the Casimir energy is specified by [4, 11, 12] 

12 

4s = 7^6^ - N v - 2) {C(5) - C(5)} , (16) 
where we have defined the function Cw(n) by 



CO 



x , cos(27r/ca;) ,„ „. 

fc=i 

Since Co,(n) < Co( n ) = C( ra ); this contribution is negative as long as N v + 2 > N H . We can 
now balance matter and gauge multiplets in exactly the same way as in the non-SUSY 
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case to ensure = 0(1). For the 2-loop contribution of the vectormultiplet, we use the 
results of Appendix B but evaluate the corresponding gaugino integrals with the shifted 
masses to obtain 

4s vector = j^yN v T(a) {((3) - C(3)} 2 . (18) 
Likewise, the hypermultiplet contribution reads 6 

<§ h V p er = -7^N v T(r) (C(3) - C(3)} 2 . (19) 



(27T) 



Since, as in the non-SUSY case, the vector multiplet contribution is positive, our 2-loop 
stabilization mechanism remains effective in this simple SUSY model as long as there is 
not too much matter charged under the most strongly coupled gauge group factor. 



4 Compactifications with branes and fixed-points 

In this section, we focus on the 5d gauge theory case discussed in Sect. El since it is 
presumably more likely to be part of phenomenologically interesting theories than the 
scalar toy model of Sect. El However, most of our discussion applies, qualitatively, also 
to the scalar case and presumably to many other 5d models with dimensionful couplings 
and corresponding 2-loop Casimir stabilization. 

As already mentioned above, the pure S 1 case discussed up to now can not give 
rise to realistic models since the 4d particle spectrum is necessarily vector-like. A simple 
way of embedding this type of S 1 stabilization in a realistic construction is to add a 
3-brane on which matter fields can live. 7 To apply our stabilization mechanism without 
any modification, one could simply assume that the brane fields are not charged under 
the bulk gauge theory responsible for the stabilizing 2-loop Casimir energy. 

If one does not make this assumption, the charged brane fields generate, at the 1-loop 
level, brane localized gauge-kinetic terms for the bulk gauge theory. Such brane-kinetic 
terms contribute to the 1-loop Casimir energy induced by bulk fields. We can view this as 
a 2-loop effect including both brane and bulk fields. It is easy to convince oneself that this 
contribution is not parametrically suppressed relative to the bulk 2-loop effect calculated 
in Sect. 01 Thus, we are precisely in the situation of 1-loop Casimir stabilization with 
brane kinetic terms considered in [3]. Our 2-loop calculation is then simply a finite 0(1) 
correction to this stabilization mechanism. However, since the coefficient of the relevant 

fi (i) (i) 

b The relation c g g vector + c g g hyper = for an adjoint hypermultiplet representation reflects the fact 
these J\f = 2 multiplets combine into Af = 4 vector multiplets. The Scherk-Schwarz mechanism leaves 
TV = 2 SUSY unbroken and the Casimir energy remains zero. 

7 As already mentioned earlier, a brane localized cosmological constant can be used to tune the 
effective 4d cosmological constant to zero. Its interplay with a possible bulk cosmological constant can 
also affect stabilization (for more options along these lines see, e.g., [13]). Here, we require the bulk 
cosmological constant to be small because of some symmetry principle (e.g. bulk SUSY) and treat the 
sum of brane cosmological constants as an unknown parameter to be fixed by the requirement of 4d 
flatness. 
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brane-kinetic term is UV-sensitive, its value after renormalization is in essence just a 
new parameter of the model. The precise result of the 2-loop calculation is then only 
meaningful if one has some first-principles knowledge about the coefficient of the brane 
kinetic term. This will, in general, require a UV completion as it is provided, e.g., by a 
string orbifold model. 

A very similar situation arises in an S 1 jZi geometry. Here, not only brane- localized 
charged fields but also the bulk gauge and matter fields themselves induce brane- localized 
gauge-kinetic terms. This is, in fact, a familiar and well-studied effect in the context of 
higher-dimensional grand unified theories (GUTs) with symmetry breaking by brane- 
localized Higgs fields and in orbifold GUTs [14,15]. It corresponds to the statement that 
a (modified) logarithmic running of the gauge couplings continues above the compactifi- 
cation scale, which can be understood as the running of the coefficient of a brane localized 
F 2 term [15]. 

In fact, it is easy to check that, even in the simple scalar model, the 2-loop vacuum 
energy in an S 1 / Z 2 compactification has (in contrast to the pure S 1 case) a logarithmic 
divergence that can be absorbed into a brane-localized 4>d 2 <f) operator. This is analogous 
to the boundary term F^ V F^ U (with //, v 6 {0, 1, 2, 3}) induced in gauge theories with 
even boundary conditions for A^. However, while in the gauge theory case only the 4d- 
part of the gauge-kinetic term is corrected (F 5fJi being zero at the boundary), in the 
scalar case the correction only affects the d§ part (the 1-loop self energy diagram being 
momentum-independent) . 

The apparent loss of predictivity associated with the UV sensitive coefficients of 
brane-localized kinetic terms is not as severe as one might naively think. The reason is 
the logarithmic enhancement of such terms associated with their logarithmic divergence. 
Indeed, the bulk gauge theory has a strong interaction scale associated with the 5d gauge 
coupling, M ~ 24ir 3 /g 2 [9]. It is natural to assume that, at this scale, the brane-localized 
F 2 term has an 0(1) coefficient. Running down to the compactification scale M c = 1/R, 
one obtains a log-enhanced coefficient ~ ln(M/M c ) = ln(247r 3 i?/ g 2 ). This is dominant 
with respect to the unknown 0(1) initial value. The corresponding log-enhanced 1-loop 
Casimir effect contribution of the brane operator is also dominant with respect to the 
true 2-loop effect calculated in Sect. El Thus, the leading piece of the g 2 /R 5 contribution 
is calculable on the basis of the low-energy field content of the model. The parametric 
behaviour is, in fact, not a pure power of R but includes the logarithmic R dependence 
of the coefficient. The 4d vacuum energy of a 5d gauge theory compactified on S l /Z 2 
thus has the form 8 



8 As already mentioned in the Introduction, higher-loop radius stabilization was also discussed in [4], 
but in a very different approach. We understand that Ref. [4] treats the 4d coupling as fundamental and 
i?-independent, resulting in loop corrections ~ \n(R)/R 4 . By contrast, we consider the 5d coupling as 
fundamental, which implies a 4d coupling ~ 1/R and hence a vacuum energy loop correction ~ ln(i?) / R 5 . 

Furthermore, a potential similar to Eq. i|20[l was derived in [6] for a scalar model with Yukawa couplings 
to brane fields. Assuming an exact cancellation of the 1-loop contributions from scalars and fermions, 
stabilization was achieved by a balancing of the 1/R 5 and ln(i?)/i? 5 terms of the 2-loop result. 




(20) 
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and can, as explained earlier, give rise to quantitatively controlled radius stabilization 
at R ~ g 2 N 3 in the case of an SU(iV) gauge theory. 

To be more specific, focus on 5d SU(iV) gauge theory on S 1 jZ% (where the gauge 
group is not broken by the orbifolding) with Nf uncharged bulk fermions and with 
charged brane fermions in a representation r at the x 5 = boundary. The logarithmic 
running above the compactification scale induces an effective brane-localized F 2 term 

C D "g^ tr ( V") ln ( Mi? ) W**) ( 8T ( r ) - f T (*)) + 5 ^ ~ * R ) (-f T («))] (21) 

at scale M c . The calculations leading to this formula are explained in detail in Sect. 3 
of [16] and can also be extracted from the earlier references [14, 15]. Here the strong 
interaction scale is M ~ 24tt 3 / (g 2 N), accounting for the large- N enhancement of higher 
loops in SU(iV) gauge theory. As far as the 1-loop Casimir energy calculation using a 
summation of KK modes is concerned, the effect of this contribution can be summarized 
by an extra contribution to the momentum-space effective action, 

^ (k 2 V - k^K) - + 2b^j (k\ v - k^k u ) , (22) 

to be used for even modes (i.e. the modes) only. Here b is defined by minus the sum 
of the two coefficients of the brane-localized F 2 terms as given in Eq. (J2T|) and k is the 
euclidean 4-momentum. The extra factor of 2 arises since the even higher KK modes are 
twice more sensitive to brane operators than the zero mode. 

Based on this, the 1-loop contribution to the 4d potential induced by the brane- 
kinetic terms reads 

V**»{R) = {ln[(l + &')^ 2 + (nM c ) 2 ] - ln[A; 2 + (nM c )]} 

3 dim(G) r d d ^k b'k 2 
~ 0. 2s I (2 7 r) d - 1 k 2 + (nM c V ' U ° j 



n=l 



Here V = 2b/(irR/2g 2 ) and the summation is only over even non-zero (cosine) modes of 
the KK mode expansion on the S 1 compact space. The prefactor of 3 accounts for the 
3 on-shell degrees of freedom of a massive vector. Note that we have simply interpreted 
the boundary operator as a correction to the energy of each separate leading-order KK 
mode. (This shifted KK mode spectrum has been derived in a different way in Appendix 
B of [17].) To be more precise, one would have to re-diagonalize the quadratic-order 
Hamiltonian after inclusion of the boundary terms and then work with the modified KK 
mode expansion [3]. However, as long as \n(MR) <C MR, we are entitled to linearize in 
the boundary term and our simplified treatment is sufficient. 

Using the calculational techniques of Appendix A, the above formula is evaluated to 
give the brane-induced correction (equivalent to the dominant part of the 2-loop Casimir 
energy) 

^w ^y , (24) 

(2tt rR °7T 
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leading to the final result 9 



V(R) 



3C(5) (4N f -3(N 2 -l)-5 iV 2 -l 
(27r) 6 i? 4 \ 2 + "(2^) F 



8T(r) 



23 



T(a) 



# 2 ln(Mi?) | 



1? 



(25) 

It is now easy to arrange for the coefficient of the l/R 4 piece to be negative and 0(1) 
while keeping the 1/R 5 term positive and 0(g 2 N 3 ) such that, as before, controlled 2- 
loop stabilization at a moderately large radius is achieved. 10 Compared with the S 1 case, 
where MR ~ iV 2 we now have MR/ \n(MR) ~ N 2 . Using asymptotic properties of the 
Lambert W function (see, e.g., [18]) the solution for N large can be written as 



MR ~ N 2 v{l + 0(\nv/v)}, v — lniV 2 
thus giving an additional enhancement factor lniV 2 . 



(26) 



5 Conclusions 

We have presented a mechanism stabilizing the size R of an extra dimension compactified 
on S 1 or on the orbifold S 1 /^ which is based on the presence of dimensionful couplings 
- a generic feature of field theoretic models with d > 4. We have shown that, by balanc- 
ing the 1- and 2-loop contributions to the Casimir energy, a perturbatively controlled 
minimum at moderately large values of R can be realized. 

In the case of scalar massless 4 theory on S 1 , the Casimir energy is calculable and 
finite as long as the masslessness is enforced as a renormalization condition. The 1-loop 
contribution scales like ~ —R~ A , while the 2-loop effect of the scalar self- interaction 
gives ~ +XR~ 5 , thus producing a nontrivial minimum. In order to ensure that the result 
is perturbatively controlled, we have added weakly coupled fermions. Their effect is 
to reduce the numerical factor of the 1-loop term while only very mildly affecting the 
higher-loop contributions. This shifts the minimum to larger values of R. 

The situation is basically the same in the case of a gauge theory. As before, the 
purely bosonic theory already produces a 2-loop minimum in the radion potential, while 
the inclusion of fermions uncharged under the most strongly coupled gauge group factor 
ensures the perturbativity of the result. We have also shown that in a SUSY extension of 
this scenario (with SUSY-breaking a la Scherk-Schwarz), the mechanism works without 
modification. 

The above calculable two-loop stabilization is modified but not destroyed if, instead 
of the circle, the orbifold S 1 /^ is considered. The crucial new feature of this situation 
are brane-kinetic terms for the gauge fields, which are generated in the 1-loop effective 

9 The factor of \ in the one-loop contribution is due to the orbifolding, since half of the modes are 
projected away (cf. also the summation in Eq. 

10 To turn over the sign in the second term, we need to introduce a certain amount of brane fermions. 
Note also that assigning negative parities to some gauge fields - in other words, breaking the gauge 
group by orbifolding - would reduce the contribution of the gauge fields to the brane kinetic terms and 
can even flip its sign. 



11 



action. Without the knowledge of the underlying UV physics one cannot predict the 
corresponding counterterm at the cutoff scale M. However, in calculating the potential 
of R for large values of the radius, R ^> M" 1 , one effectively integrates out the physics 
down to that scale, and the logarithmic running of the brane-kinetic term dominates 
over its unknown initial value. The dominant 2-loop effect can then be obtained by 
evaluating the 1-loop integral in the presence of this log-enhanced brane-kinetic term. 
The new contribution scales like ~ g 2 ln(R M) / R 5 and is calculable on the basis of the 
bulk and brane field content of the model. We emphasize that this log-divergent part 
of the 2-loop calculation dominates over the remaining (finite) 2-loop effects (scaling as 
~ 1/-R 5 ) as well as the one- loop effect of the unknown brane-kinetic term (whose leading 
effect for large R is also ~ 1/-R 5 )- 

Although we have focussed exclusively on 5d models, we note that our mechanism is 
equally suitable for other higher-dimensional manifolds with a single unstabilized mod- 
ulus. This is because such theories look very similar from a 4D point of view, the main 
difference being a modified KK-mass spectrum which depends on the value of the modu- 
lus. Various possibilities for stabilizing such a single modulus have recently been discussed 
in the context of the KKLT proposal [19], and we believe that Casimir energies will be 
relevant in this context. This has, in fact, very recently been discussed in the context 
of the 1-loop effect of massive vector fields in [20], and one can expect that the higher- 
loop contributions discussed here will also play an important role in the further study of 
moduli stabilization in more complete, higher- dimensional scenarios. 

As a simple and specific example, we briefly consider S n compactifications with the 
volume undetermined by the Einstein equations. The Casimir energy coming from a self- 
interacting scalar field will take the same general form as in Eq. ©, with the powers of 
R modified according to the dimensionality of the coupling and the cM to be calculated 
from the appropriate KK spectrum. The coefficient c^ 1 ' has been calculated for the case of 
a sphere (without additional infinite dimensions) in Ref. [21] for n— 1,2,3, 4 and found 
to be negative. Moreover, it is clear that C ( 2 ) is always positive since the integral appears 
squared. Thus, if the inequality c^ 1 ' < survives the transition from S n to M4 x S n 
for n > 1 (as it does for n = 1), our stabilization mechanism extends straightforwardly 
to these geometries. In any case, we can be confident that many examples of 2-loop 
(or higher-loop) Casimir stabilization exist within the rich class of models with n > 1 
compact dimensions where the total volume is a modulus at tree level. 

Our final point concerns SUSY theories on S l /7*2- The prototype SUSY breaking 
mechanism on S 1 , the Scherk-Schwarz mechanism [10], extends to this case. It can equiv- 
alently be described by giving a vacuum expectation value to the F component of the 
radion superfield (the chiral superfield whose lowest component contains R) [22,23]. 
Indeed, the resulting tree level action corresponds to no-scale supergravity and the po- 
tential for R is completely flat. Stabilization at 1-loop has previously been studied by use 
of bulk mass terms for hypermultiplets [3,12,24], brane kinetic terms [3,25] or massive 
vector multiplets [20]. The mechanism we described in Sec. 0] seems to be quite suitable 
for SUSY S ,1 /Z 2 models without 5d masses. The 1-loop contribution is given by ^ of the 
S 1 result obtained in this paper. We expect the dominant 2-loop effect to be the log- 
enhanced brane-kinetic terms, the same as in the non-SUSY case. However, a detailed 
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analysis of this scenario, of its interplay with the other stabilization mechanisms men- 
tioned above, or even the construction of a realistic SUSY model go beyond the scope of 
the present paper and are left to future research. 
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Appendix A: Basic 1- and 2-loop integrals 

We begin by rederiving the known 1-loop vacuum energy for a real scalar field on R A x S 1 
in dimensional regularization. The 1-loop effective potential of the 5d euclidean theory 
is given by the d — > 5 limit of 

.. +00 „ 7(2—1 7 

V PM = -jR E TJiHk 2 +(nM c f}, (A.l) 

n=— 00 ^ ' 

where M c = 1/R. Using the fact that the ci-dimensional integrals of any power of k 2 and 
of ln(/c 2 ) vanish in dimensional regularization, we have 



M c f M - 2 fd d - l k k 2 7rcoth(7r|A;|/M) 







iM *Wfm \k\/M ■ < A - 3 > 



Appealing again to the fact that the d- dimensional integral of a pure power vanishes, 
coth(7r|/c|/M) can be replaced by coth(7r|A;|/M) — 1, making the d d ~ 1 k integral finite and 
allowing one to take the limit d — > 5. The resulting 4d effective potential is 

V£\R) = 2kRV^\r) = -J®, with C(5) = 1.0369... , (A.4) 

in agreement with Eqs. Q and (jlOj). 

Higher-loop corrections to the effective action (in particular to — V^) are given by 
the sum of all one-particle-irreducible vacuum diagrams of the relevant field theory. As 
already pointed out in the main text, the 2-loop correction is essentially just the "figure 
8" diagram. To be more precise, one has to add to the lagrangian the mass counterterm 
of the uncompactified 5d theory ensuring that the scalar remains massless at the 1- 
loop level. This counterterm becomes part of a tadpole-like 1-loop vacuum diagram, 
which is part of the 2-loop correction. However, in dimensional regularization the mass 
counterterm vanishes and the 2-loop correction (for a single scalar field) simply reads 

T,(2) A 2 _ 1 ^ f d d ~ l k 1 E , 

V* = 77 1 , where / = — > / - — -7 — . -— . (A.5) 

5 8 R f-^ J {2ix) d k 2 + (nM c ) 2 K ' 
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Here the prefactor contains a symmetry factor (1/8) and the tadpole integral I is eval- 
uated as the d —>■ 5 limit of 

f d^k TrcothMkl/M) C(3) . , , nrtn , , A „* 

This completes the derivation of Eq. (JTTJl and thereby of Eq. (JTJj). 



Appendix B: SUSY checks of gauge theory results 

To verify the results of our 2-loop gauge theory calculation, we have performed checks 
based on the vanishing of the Casimir energy in theories with unbroken SUSY in 5 and 
4 dimensions. 

Specifically, we can consider a 5d super Yang-Mills (SYM) theory which, in addition 
to the gauge fields, includes a real adjoint scalar E and a Dirac gaugino. To calculate the 
2-loop effect in this theory, coupled to a hypermultiplet (a Dirac fermion and two complex 
scalars) in a representation r of the gauge group, we need the contributions coming from 
the additional Yukawa- and scalar self-interaction diagrams shown in Fig. EH For the 
explicit lagrangian see, e.g., [26]. 



Figure 3: Two-loop diagrams from Yukawa couplings and scalar self-interactions. 
The three "setting-sun" diagrams coming from the Yukawa couplings yield 

J 2 ) _ _o (2) J./ \ 

gaugino— fermion— scalar V / > 

4 2 L,n- E = -4 2) T(r), (B.l) 
c (2) - -c {2) T(a) 

c gaugino-S ~~ °0 1 \ a ) • 

Furthermore, in the scalar sector there are two new "figure-8" diagrams coming from the 
coupling of £ to the scalar and of the scalar to itself, giving two new contributions 

C scalar-£ = ^ C V ^X r ) > 
C Lalar-scalar = ^ Cq ^ T(r) . (B.2) 

Finally, the gaugino and S as well as the hypermultiplet fermion and scalar are charged 
under the gauge groups. The corresponding 2-loop corrections can be directly read off 
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from Eq. (|15|). The two bosonic contributions are 

c (2) - -T(a)c {2) 

L vector-E — ^ J \ a ) c ) 
C vector-scalar = ^( r ) C o'i (B-3) 

where we have included factors \ and 2 to account for the reality of E and for the presence 
of two hypermultiplet scalars respectively. The two fermionic contributions are precisely 
as in the Eq. (|15|). just with T(r) replaced by T(a) in the case of the gaugino. 

The contributions of the charged hypermultiplet thus sums up to 

c (2) +c {2) +c (2) +c {2) +c (2) +c {2) 

vector— scalar 1 scalar— £ 1 scalar— scalar 1 vector— fermion 1 fermion— £ 1 gaugino— fermion— scalar 

= (7 + 2 + 3-3- l-8)T(r)4 2) = 0. (B.4) 
Similarly, the self-interactions of the vector-multiplet give rise to 

c (2) + c (2) + c {2) + c (2) -(- + --?, -l\ T(a) c {2) - (B 5) 

^vector i C vector-S ' C vcctor-gaugino ' L gaugino-£ 1 4 4 / ^ ' _ l" ' / 

Next, since we have kept the d dependence in Eq. f)15|) . we can immediately extend 
our analysis to a 4d M = 1 SYM theory. The 24oop Casimir energy contribution is 
proportional to 

C 4d vector + C 4d vector-gaugino = — l)^( a ) C ^ = ^ ' (B-6) 

where the fermionic contribution includes an extra factor | relative to Eq. (fTK|) to account 
for the chirality of the 4d gaugino. 

Finally, consider the contribution of charged 4d matter. The 24oop effects based 
on gauge interactions can be inferred from Eq. (|15j) (with an appropriate factor | for 
fermion chirality in the loop). The gaugino- fermion-scalar contribution is as in Eq. (|B.1|) . 
but with an extra factor | for the fermion chirality and half the number of scalars. The 
scalar "figure 8" diagrams are induced by the D term potential and give rise to | of the 
contribution displayed in Eq. (jB.2|) because of the missing trace over SU(2)r degrees of 
freedom, XxoiOi = 6. 11 Overall, one finds 

c (2) +c (2) +c (2) +c (2) 

4d vector— scalar 1 4d scalar— scalar 1 4d vector— fermion 1 4d gaugino— fermion— scalar 



5 1 
2 + 2 



>W)c< 2 > = 0, (B.7) 



concluding our SUSY based checks. 



11 The scalar "figure-8" diagram also contains a piece proportional to (trT a ) 2 which is zero for non- 
abelian as well as anomaly-free abelian gauge theories. The reason why the Casimir energy does not 
vanish in the presence of an anomalous U(l) is the associated occurrence of a one-loop Fayet-Iliopoulos 
term which spontaneously breaks SUSY. 
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